PCCF-RI-0414 
ECT*-04-10 



Direct CP Violation in B decays with p — u Mixing 

Ziad J. Ajaltouni 1 
Laboratoire de Physique Corpusculairc 
Universite Blaise Pascal CNRS/IN2P3 
F63177 AUBIERE CEDEX FRANCE 

Abstract 

A complete study of the processes B — » n + iv~V (V = 1 ) is performed both in the framework 
of the helicity formalism and the effective lagrangian approach. Emphasis is put on the factorization 
hypothesis and the importance of the p° — uj mixing in enhancing the direct CP violation. New results 
involving some branching ratios and the ratio of the Penguin/Tree amplitude are given in details. 

1 Physical Motivations for the channels B — > V1V2 and their De- 
cay Kinematics 

In the framework of the LHCb experiment devoted to the search for CP violation and rare B decays, 
special care is given to the B decays into two vector mesons, B — ► V1V2, V%= 1 ■ 

(i) B decays being governed by weak interactions, the vector- mesons are 'polarized and their final states 
have specific angular distributions; which allows one to cross-check the Standard Model (SM) predictions 
and to perform tests of models beyond the SM. 

(ii) In the special case of two neutral vector mesons with V — V° ; orbital angular momentum t, total 
spin S and CP eigenvalues are related by the following relations : t = S = 0, 1, 2 CP = (— l) e 
, which implies a mixing of different CP eigenstates, proving a CP non-conservation process. According 
to Dunietz et al £Q, tests of CP violation in a model independent way can be performed and severe 
constraints on models beyond the SM can be set. Because the B meson has spin 0, the final two vector 
mesons, V\ and V 2 , have the same helicity Ai = A2 = — 1, 0, +1, and their angular distribution is isotropic 
in the B rest frame. Let H w be the weak Hamiltonian describing the B decays. Any transition amplitude 
between the initial and final states will have the following form: 

H x = (V^V^H^B) (1) 

where the common helicity is A = — 1,0, +1. Then, each vector meson V% will decay into two pseudo- 
scalar mesons, at and bi ; ai(bi) can be either a pion or a kaon which angular distributions depend on Vi 
polarization. 

The helicity frame of a vector-meson Vi is defined in the B rest frame such that the direction of the 
Z-axis is given by its momentum pi . Schematically, the whole process gets the form: 

B — ► Vi + V 2 — ► (ai + 61) + (02 + b 2 ) ■ 

The corresponding decay amplitude, M\(B — > X)i=i( a « + ^i))> ^ s factorized according to the relation, 

2 2 
M X (B -> 5^(04 + h)) - H X (B -> Vi + Vt) x YlMV -► a, + b. t ) , (2) 
i=i i=i 

where the amplitudes Ai(Vi — > a, + bi) are related to the decay of the resonances Vi. For a given value of 
A and a well defined final state, amplitudes Ai(Vi — > Oi + bi) are given, according to the Wigner-Eckart 
theorem, by the following expressions: 

A x (Vi ->oi + &i) =c 1 D 1 Xmi (0,6 1 ,0) and A 2 (V 2 -> a 2 + b 2 ) = c 2 D\ m2 (0, 2 , -0) . (3) 
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In Eq. (J3J), the c\ and C2 coefficients represent respectively the dynamical decay parameters of the 
V\ and V2 resonances. The term D\ {(j)i,6i,—(f>i) is the Wigner rotation matrix element for a spin-I 
particle and we let A(ai) and \(bi) be the respective helicities of the final particles ai and bi in the Vi rest 
frame. 8\ is the polar angle of a\ in the V\ helicity frame. The decay plane of V\ is identified with the 
(X-Z) plane and consequently the azimuthal angle 4>\ is set to 0. Similarly 82 and <p are respectively the 
polar and azimuthal angles of particle 02 in the V2 helicity frame. Finally, the coefficients m,i are defined 
as mi = A(a,) - \{bi). 



2 Decay Dynamics and Basis for Simulations 

The most general form of the decay amplitude A4(B — > Y^i=i( a i + bi)) 1S a ^ near superposition of the 
previous amplitudes M\(B — > J2i=i( a i + M) denoted by: 

2 2 

i=l A i=l 

The decay width T(B — ► V1V2) can be computed by taking the square of the modulus, \M.(B — ► 
Si=i( a » w hich involves the three kinematic parameters, 81,62 and 0. This leads to the following 

general expression: 

1 2 2 

d 3 r(5 ^ Vi^) cx 5^M A (B-f5^(oi + fci)) =^/u,A'F A , A ,(0 1 )G A , A ,(0 2 ,0) , (5) 

A i=l A, A' 

which gives rise to three density-matrices h\\> , F AA ' (#i ) and G AA ' (#2 , 4>) '■ (i) The factor h\\> = H\Hy 
is an element of the density-matrix related to the B decay; (ii) -F AA ' (8\ ) represents the density-matrix of 
the decay Vi — ► ai + b\ and (iii) G\y(02, 4>) represents the density- matrix of the decay V2 — » &2 + 62. 
The analytic expression in Eq. JSJ exhibits a very general form. It depends on neither the specific 
nature of the intermediate resonances nor their decay modes (except for the spin of the final particles). 

The previous calculations are illustrated by the reaction B° — > K*°p° where K*° — > K + tt~ and 
p° — ► 7r + 7r~. In this channel, since all the final particles have spin zero, the coefficients mi and m2, 
defined previously, are equal to zero. The three-fold differential width has the following form : 

d 3 T{B ViV 2 ) _ + h __j sin 2 0lSin 2 62 / 4 + h^cos^cos^ 



d(cos 9i)d(cos 8 2 )d(j) 

+ |^(/i+o)cos <ft — J 2r m(/i+o)sin cf> + 3$e(ho-)cos <j> — J%i(/io-)sin cf> |sin 2#isin 202/4 

cos 2</> — ^fm{h-\ )sin 2</)|sin 2 0isin 2 02/2 . (6) 

It is worth noticing that the expression in Eq. JBJ is completely symmetric in 8\ and 82 and consequently 
the angular distribution of a\ in the V\ frame is identical to that of 02 in the V2 frame. From Eq. JSJ the 
normalized probability distribution functions (pdf) of 8\, 82 and (/) can be derived and one finds : 

/(cos0i, 2 ) = (3h o - I)cos 2 6q !2 + (f - h 00 ) , g{<j>) = I + 2 &e(h+- )cos 2<p - 2 J%i(/i+_ )sin 2<p . (7) 

A practical way to compute the matrix elements is to use the Effective Hamiltonian approach based on 
the general hamiltonian : 

H eff = 7§ E VckmC^O^) , (8) 

where Gf is the Fermi constant, Vckm is the CKM matrix element, Cj(/x) are the Wilson Coefficients 
(W.C.), Oi(/i) are the operators associated to the tree, QCD-penguin and EW-penguin diagrams 
and, finally /1 is the renormalization energy scale taken equal to mg. 
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Then, applying the Operator Product Expansion (OPE) method pioneered by Wilson, the W.C. Cj 
are calculated perturbatively at the Next to Leading Order (NLO) for an energy scale > tub [IJ- The 
non-perturbative effects which are related to the operators Oi and representing physical processes at an 
energy < ms are introduced through a set of form factors. The latters are explicitly computed in the 
framework of the pioneering BSW models However some free parameters remain like : (i) the ratio 
q 2 /nib 2 where q 2 is the squared invariant mass of the gluon appearing in the penguin diagrams and (ii) 
the effective number of colors N c e ^ . 

Final State Interactions and p° — u> Mixing 

Hadrons produced from B decays are scattered again by their mutual strong interactions, which could 
modify completely their final wave- function. Computations of the branching ratios B — > hadrons must 
take account of the final state interactions (FSI) |3] which are generally divided into two regimes : 
perturbative and non-perturbative. These two aspects have been already mentioned above in the framework 
of the OPE method. However an important question arises: how to deal with the FSI in a simple and 
practical way in order to perform realistic and rigorous simulations? 

The method which has been followed for the computations is largely developed in |2] and |S] and it is 
based on the hypothesis of Naive Factorization, which can be summarized as follows : 

★ In the Feynman diagrams describing the B decays into hadrons like tree or penguin diagrams, the soft 
gluons exchanged among the quark lines are neglected. 

★ The color number N c is no longer fixed and equal to 3. It is modified according to the relation : 
( N *ffj = f + £ where £ is an operator representing the non-perturbative effects. 

★ The QCD-penguin diagram introduces an intrinsic phase-shift, Sp/? , by comparison with the tree one 
(BSS mechanism [7]). Thus, the total amplitude gets an absorptive part, which is an illustration of the 
FSI in the perturbative regime. 

★ Another important effect which appears in the channels B — > ir + Tr~V is the p° — u> mixing, which is 
an unavoidable quantum process. Indeed, the tree amplitude A T and the penguin one, A , are modified 
according to the following relations : 

(K*n-ir+\H T \B) = -^-fi^A, + , (K*tt-it+\H p \B) = -^-fi^pk, + % p (9) 

SpSoj Sp SpSijj Sp 

Here ty (V — p or lu) and py are respectively the tree and penguin amplitudes for producing a vector 
meson V , g p is the coupling for p° — > ir + ir~ , tl puj is the effective p° — u> mixing amplitude and sy is the 
inverse propagator of the vector meson V , sy = s — m,y + imyTy where y/s is the invariant mass of 
the 7r + 7r~ pair. 

The ratio A p / A T , which is a complex number, gets the final expression : 

re'V* = n p" p " + SqjP p ) (10) 

Hpcj^cj ~t~ Sojtp 

where 5 is the total strong phase arising both from the p° — ui resonance mixing and the penguin diagram 
quark loop, and <j> is the weak angle resulting from the CKM matrix elements. 

3 Main Results and Comparison with Recent Experimental Data 

Owing to the presence of resonances with large widthes, the mass of each resonance is generated according 
to a relativistic Breit- Wigner distribution : 

da T R M R 

oc 



dM 2 (M 2 - M 2 f + (T R M R ) 2 ' 
Mr and Tr being respectively the mass and the width of the vector meson. 

Then, combining both the Wilson Coefficients and the BSW formalism and including the p° — u mixing, 
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the helicity amplitude (computed in the B meson rest-frame) is given by the following expression : 
H X (B - p°(u)V 2 ) = iB p x (V ub V* s c p tl - V tb V t * s <f p J + *C7£(K*V£< - V tb V*c? 2 ) + 

■p \iB%(V ub V: s cl ~ VaVZqj + iCt{V ub V: s cl - V tb V t * s c% 2 )} , (11) 



(s p - ml) + in\ 

where the terms B^ and C x are combinations of different form factors. Their explicit expressions, 
corresponding to the helicity values (A = — 1, 0, +1), are given in Ref. This expression allows to 
deduce the dynamical density-matrix elements h\y given by : 

^AA' =H X (B^ P °(lu)V 2 )H* x ,(B -> p°(u)V 2 ) . 

Because of the hermiticity of the DM, only six elements need to be calculated. The main results are : 

1) The matrix elements hij depend essentially on the masses of the resonances. Their spectrum of is 
too wide because of the resonance widths, especially the p° width T p = 150 MeV/c 2 . 

2) The longitudinal polarization, h 00 = \H Q \ 2 , is largely dominant. In the case of B° — > p°(u})K*°, 
the mean value of h 00 is as 87% while for B + — > p°(u;)p + , its mean value is as 90% (Fig.l). These results 
have been confirmed recently by both BaBar [Sj and Belle collaborations UJ. 

3) The matrix element = |-ff-i| is very tiny, < 0.5%. 

4) The non-diagonal matrix elements hij are mainly characterized by : (i) The smallness of both 
their real and imaginary parts, (ii) The ratio ^m/^e as 0.001 — * 0.1. (iii) In the special case of 
B + — > p°(uj)p + , c^m(hij) as 0.0 . Our conclusion is that there is a kind of universal behavior of the 
density-matrix elements, whatever the decay B — > tt + tt~ -V is (V = K*°, K*^ , p±). 

• Consequences on the Angular Distributions : 

In the helicity frame of each vector-meson Vi , the angular distributions given above (see Eq. (7)) become 

simplified : because of the small value of ), the azimuthal angle distribution g((f>) is rather flat. In 

the expression of /(cos^i^) , the longitudinal part hoo being largely dominant, the polar angle distribution 
is as cos 2 #. 

Branching Ratios and Asymmetries 

★ The energy Ei and the momentum pi of each vector meson vary significantly according to the generated 
event. So, the width of each channel is computed by Monte-Carlo methods from the fundamental relation : 



dT(B — > ViV 2 ) - ^jjlMiB^VM)] 2 ^^±5\P- Pl -p 2 ) (12) 



1tt 2 M 1 { 1 2I] 2E 1 2E 2 

from which the specific branching ratios are deduced. 

* For a fixed value of q 2 /m b 2 , the BRs depend strongly on the Form Factor models. They could vary up 
to a factor 2. 

* The relative difference between two conjugate branching ratios, Br(B — » /) and Br(B —>■/), is almost 
independent of the form-factor models. The global asymmetry defined by : 

_ Br(B -> f) - Br(B -> /) 
CP Br(B -» /) + Br(B -» /) 

is usually < 2% 

* However, an interesting effect is found in the variation of the differential asymmetry with respect 
to the 7T7T invariant mass. This parameter defined as : 

aCP(m) = T m (B^f) + r m (B->f) 

is amplified in the in the vicinity of the u> resonance mass (±20 MeV/c 2 around M w = 782MeV/c 2 ). 
acp{rn) is as 15% in the case of B° — * K*°p (uj) and equal to 24% in the channel B ± — * p ± p°(uj). 

This kind of asymmetry is almost independent of the form factor models. It is worth noticing that 
this novel effect has been predicted analytically in the channel B — > VP — * ir + ir~ir by Leitner et al 
and its only explanation is the mixing process of the two vector- mesons p° and u. 
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Ratio Penguin/Tree 

The ratio Penguin/Tree is given by the following relation derived from equation (11) : 

P iS icb n 'tbVtq i 

T= re e ' r = r \y^ 1 

J v ubV uq 

where r' is the " naked" ratio P/T . It is almost constant over the tttt interval mass, but it varies very 
sharply in the ui interval, from 760 MeV — > 820 MeV, especially in the channel B° — > K*°p°(ui) where 
it reaches 60%. Its variation is almost independent of q 2 /m\ but it depends on N%ff. 



Channel 


Usual Values 


a; Interval 


B° -> K*°p°(uj) 


0.08 < r' < 0.30 


0.60 


B+ -> K* + p u (w) 


0.04 < r' < 0.05 


0.06 


B+ -> pV'M 


r' w 0.016 


0.05 



Final Phase-Shift J 

The strong phase 5 which is the phase difference between the Penguin and Tree diagrams is the main 
ingredient of the absorptive part of the B decay amplitude. Its physical origin is related to the : (i) Intrinsic 
pahse-shift induced by the Top quark in the Penguin diagram, (ii) the complex Wilson Coefficients, and 
(hi) essentially the p° — u> mixing in the tttt final state interactions. It depends on the ratio q 2 /rrib 2 
and strongly on N**1 ' . Usually, 5 is almost constant in all the tttt invariant mass interval except in the 
uj resonance window, 770 -> 790 MeV, where it undergoes a variation of 80° -> 100° in the K*°p°(uj) 
channel and a variation of 5° — > 25° in the K* + p°(uj) one. 

• Very interesting physical consequences can be inferred from the exhaustive study of the parameters 
P/T and 5 with the tttt invariant mass. The direct CP asymmetry parameter is defined according to the 
relation : 

dir = A 2 - A 2 = -2 r sin 5 sin^ 

CP A 2 +A 2 l + r 2 +2rcos<5cos0 ' 1 1 

where $ is one of the weak mixing angle deduced from the CKM matrix elements. In the channel 
B — > p°{u)K* , angle $ is identified with Arg[VtbV t */V u bV* s ] = 7 ; while in the channel B — > p°(us)p, 
angle $ is given by Arg[V t bV t * d /V u bV* d ] = (3 + 7 = II — a. So, the theoretical knowledge of r and 6 and 
the experimental measurements of af^jp according to the tttt invariant mass allow to extract angle(s) $ 
from the above equation (fT^I . 

These results could be seen as experimental challenges for the future LHC experiments in the held of 
B physics like LHCB. 



Recent Experimental Results 

Recently, B factories like BaBar and Belle experiments published interesting results related to the charm- 
less decays B — * V1V2. They both agree on the fact that the longitudinal part of the decay amplitude 
is very dominant, which is one of our essential results. However, these collaborations do not take into 
account the process of p° — uj mixing in the estimation of the branching ratios and the asymmetries. By 
computing the branching ratios from relation l|12|l and comparing them with those published in ref. [S] 
and HI, we can summarize the main results in the tabled 



4 Conclusion and Perspectives 

Helicity formalism has been used very successfully for a full computation and numerical simulations of 
the channels B — > tt + tt~ V with V = K*°, K* ± , p ± . Naive factorization is very useful for weak hadronic 
B decays despite its theoretical uncertainties. Furthermore, interesting results have been obtained like : 
(i) The important role of the form factor models, (ii) The longitudinal polarization is largely dominant, 
whatever the form factor model, (iii) The p° — uj mixing is the main ingredient in the enhancement of the 
direct CP violation, (iv) A new way to look for direct CP Violation is found and it can help to develop 



5 



Channel 


Br(xl0~ 6 ) 


h = \H \ Z 


A C p 


p K* + (BaB&T) 


10.6+^±2.4 


0.96+^|±0.04 


0.20t^±0.04 


Our results 


2.3 -» 5.8 


87% 


-6.4% -» -22% 


p°p + (BaBar) 


22.5+^±5.8 


0.97+o;^±0.04 


-0.19±0.23±0.03 


p°p+ (Belle) 


31.7i^±7.1 


0.95±0.02±0.11 


0.00±0.22±0.03 


Our results 


11.0 -> 20.0 


90% 


-8.5%-* -10% 



Table 1: Predicted results compared to experimental data from BaBar and Belle 



new methods for measuring the angles 7 and a. 

What remains to be done is to cross-check these predictions with experimental data coming soon from 
the LHC exoeriments. 

Acknowledgments : Z.J. A. is very indebted to the organizers of the QFTHEP04 conference which was 
held in this historical and marvellous city of Sant-Petersburg. 
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